arXiv:1509.01002v3 [cond-mat.str-el] 17 Mar 2016 


Nature of Possible Magnetic Phases in Frustrated Hyperkagome Iridate 
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Based on Kitaev-Heisenberg model with Dzyaloshinskii-Moriya (DM) interactions, we studied 
nature of possible magnetic phases in frustrated hyperkagome iridate, Na4lr308 (Na-438). Using 
Monte-Carlo simulation, we showed that the phase diagram is mostly covered by two competing 
magnetic ordered phases; Z 2 symmetry breaking (SB) phase and Ze SB phase, latter of which is 
stabilized by the classical order by disorder. These two phases are intervened by a first order phase 
transition line with Zg-like symmetry. The critical nature at the Ze SB ordering temperature is 
characterized by the 3D XY universality class, below which U(l) to Ze crossover phenomena appears; 
the Ze spin anisotropy becomes irrelevant in a length scale shorter than a crossover length A* while 
becomes relevant otherwise. A possible phenomenology of poly crystalline Na-438 is discussed based 
on this crossover phenomena. 

PACS numbers: 


I. INTRODUCTION 

An interplay between spin and orbital degree of 
freedom enriches the physics of Mott insulatorsi. A 
strong relativistic spin-orbit interaction in cubic IrOe 
octahedra endows Ir electron with spin-orbit entan¬ 
gled Kramers doublet (pseudo-spin doublet)^ which 
opens a new root to “J = Mott insulators with 
coordinate-dependent pseudo-spin-anisotropic exchange 
interactionsMott insulating materials with bond- 
dependent spin-anisotropic interactions possibly stabi¬ 
lize gapless quantum spin liquid state with exotic exci¬ 
tations such as Majorana fermions^. Honeycomb-lattice 
Na2lr03^~— a, /3, 7-Li2lr03^ii^ii^ and hyperkagome- 
lattice Na4lr308 (Na-438)i2i“— are being intensively ex¬ 
plored along this research interest. Throughout ther¬ 
modynamic measurement s^ d^d^d^ resonant inelastic X- 
and neutron scattering experiments^iii^, the 
nature of low-temperature magnetic phases in all the 
honeycomb iridate compounds has been mostly clar¬ 
ified. Meanwhile nature of a low-T phase of Na- 
438 is still veiled in mystery, although possible mag¬ 
netic nematio2ii2i, valence bond solid ^ and spin liq¬ 
uid phases22ii2^ are being proposed theoretically. Early 
thermodynamic measurements such as magnetic specific 
heat, heat capacity, and magnetic susceptibility indicate 
spin liquid feature in Na-43 8^'^i^^ , while recent neutron 
scattering and muon spin relaxation experiments on pow¬ 
der samples suggest possibilities of a short-range ordered 
quasistatic spin state (dubbed as “configurationally de¬ 
generate phases with fluctuating order”) or spin-freezing 
like phase in Na-438 One of the major obstacles 

toward comprehensive understanding of low-T magnetic 
properties of Na-438 stems from its low crystal symmetry, 
which results in complexity of electronic band structure 
and effective exchange model of localized spins 

In this paper, we introduce an effective spin exchange 
model for the hyperkagome iridate to obtain a compre¬ 
hensive understanding of classical magnetism possible in 
Na-438. Based on a lattice parameter of Na-438, we 


postulate additional lattice symmetries other than an 
exact crystal symmetry, to derive a relatively simpler 
but realistic effective spin model (Sec. II). Using Monte 
Carlo (MC) simulation and Luttinger-Tisza (LT) analy¬ 
sis, we derive a classical magnetic phase diagram for Na- 
438 (Sec. III). The phase diagram is mostly covered by 
two phases; one is Z2 symmetry breaking (SB) magnetic 
phase and the other is Ze SB magnetic phase. The Ze 
anisotropy in the latter phase is attributed to the entropy 
effect (Sec. IV). At finite temperature, these two phases 
are separated by a first order phase transition line with 
Zg symmetry. The finite size scaling (FNS) analysis con¬ 
cludes that criticality at the ordering temperatures of Z2 
and Ze phases are characterized by the 3D Ising, and 3D 
XY universality class respectively. For a finite-size sys¬ 
tem, an intermediate temperature regime appears below 
the ordering temperature of Ze phases, where the Ze spin 
anisotropy becomes effectively irrelevant and spin order¬ 
ing develops in a U(l) symmetric way. This crossover 
temperature regime AT* = Tc — T* is scaled with the 

_ i_ 

linear dimension of the system size A as AT* K 
with vq = 1.45 ^ 1.85 (Sec. V). Effects of the quantum 
fluctuation are also mentioned in Sec. VI. A possible phe¬ 
nomenology of low-T magnetic behaviors of the powder 
samples as well as a brief summary are given in Sec. VII. 


II. AN EXCHANGE SPIN MODEL FOR 
HYPERKAGOME IRIDATE 

The hyperkagome lattice is a three-dimensional lattice 
which comprises of corner-sharing triangles (Fig.[T]). The 
cubic unit cell contains 12 crystallographically distinct 
lattice points. Each sublattice point has a two-fold rota¬ 
tional axis, around which the lattice is symmetric under 
the C 2 rotation; C 2 @j {j = 1,2-•• ,12) and j sublat¬ 
tice index. We can choose C 2 @j {j = 1, • • • ,6) as gen¬ 
erators of the exact crystal symmetry group; the others 
six are identical to one of the generators {C 2 m = ^202, 
C2@8 = 672@i, C2@9 = C2@6^ 672@10 = 672@5, C2@ll = 
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<^2035 C' 2 @i 2 = <^204 in Fig. [T]). The two-fold rotational 
axes of C 2 @i^ ^ 2@2 7 ^ 2@3 7 ^ 2@4 5 ^ 2@5 nnd C72@6 nre along 
(0,1,1), (1,1,0), (1,0,1), (1,-1,0), (0,1,-1), and (1,0,-1) di¬ 
rections in the cubic unit cellj^i^. By these C 2 rotations. 


a form of the exchange interaction between a pair of near¬ 
est neighboring Ir pseudo-spin doublets determines all 
the others; 


J 


/ Ji Gi2 + Di2 G 31 + 1^31 \ / \ 

H = J2{ Si,, Si,p ) C?12 - £>12 J2 G23 + £>23 Sj,, , 

{i,j) \ G31 — Dsi G23 — D23 J3 J \ Sj^p j 


( 1 ) 


with 


( ~ 

{Siij),x ^Siij),y ^i{j),z) for (i,j) = (1,2), (10,12), 
( (bj) ~ (2,3), (4,9), 

( (bj) ~ (3,1), (b,5), 

( ^Kj)iy j (bj) ~ (1", 3), (12,1), 

( ^Si^j),z TSi(j)^y ) for (i,j) = (4, 5), (2,10), 

^ ( ^iU),y ) foi* (bi) = (9,2), (11,7), 

( (bj) ( 3 , 3 ), ( 10 , 9 ), 

( ^^i{j),y ^^i{j),x ^i(j),z ) F)r {i-fj) (3,4), (6,12), 

( ^^i{j),x ^i(j),y ^^i(j),z ) F)r {i-jj) (8,4), (5,7), 

( TSi(j)^z Siij),x ^^iij),y ) foi* (^G*) = (7,6), (12,11), 
( ^^iU),y TSi(j),x ) for (ij) = (9,8), (11,10), 

^ ( ^^i{j),x ^i{j),z ^i{j),y ) fo^ (bj) (1, 9), (8, 11), 

= I,--' ,12 the sublattice index and the summation 
of (i, jf) is taken over all the nearest neighbor Ir sites. The 
exchange spin interaction takes the 3 by 3 real-valued 
matrix form, containing both symmetric {Ji and Gij) 


and antisymmetric component (Dij). 

Quantitatively, spin anisotropy in hyperkagome iridate 
is determined by relative strength between the atomic 
spin-orbit interaction and non-cubic crystal field split¬ 
ting among t 2 g orbitalsIn the larger spin-orbit cou¬ 
pling case, Jeff = ^ doublet respects the cubic symmetry 
of the crystal field To reduce a number of the spin- 
model parameters in eq. o, we will further postulate 
the cubic symmetry of the crystal field as well as several 
additional symmetries which are nearly respected by the 
lattice parameters of Na-438. 

To this end, note first that a form of the quadratic spin 
Hamiltonian between nearest neighboring Ir pseudo-spin 
doublets (for clarity of the explanation, choose Irl and 
Ir2 in Fig. [T] henceforth) is mostly determined by direct 
transfer integral between the two doublets and indirect 
transfer integrals mediated by oxygen p orbitals neigh¬ 
boring to the two doublets (01 and 02 in Fig. [2]). An 
electronic Hamiltonian for these four (Irl, Ir2, 01 and 
02 in Fig. [2]) is given by 


— ^lr,d fj,afj,o- + 

^rn,p pln,a 

^aP'm,a,a + ^ ^ 

^2j,crcr'/i|cr/j,cr' 

ji=Irl,Ir2,cr 

m=01,02,a a=x,y,z 

i,ji=Irl,Ir2,cr,cr' 


+ E 

^im,a,aa' fi^crP'm,a,a' T H.C. 

T. fU 

],ct' 

i=Irl,Ir2,m=01,02,cr,cr' 


j=lrl,lr2,a,a' 


+ E 

Um,ab p\n,b,cFP\n,a,cF'P^^(^^^'P^J^^ 


m=01,02,a,b=x,y,z,a,a' 





+ ^ ^ Um,dp 

i=lrl,lr2,m=01,02,a=x,y,z,a,a' 


( 2 ) 


where Jeff = ^ pseudo-spin doublet respects the cubic 
symmetry^-— 

/j,± — "^73 1 ^'^^3iXZ,s 1=^ T ==F i =± I*-(3) 

^Om,p is an effective atomic energy for the oxygen Om 
(m = 1, 2) and we assume that three p-orbitals at the re¬ 


spective oxygen Om take the same atomic energy, while 
^oi,p 7^ eo 2 ,p- Sim,a,aa' stauds for the transfer between 
the Ir doublet at i = Irl or Ir2 and the neighboring oxy¬ 
gen p-orbital (a = x^y^z) at m = 01 or 02. cr(cr') is 
the pseudo-spin index, is the transfer between 

two nearest neighbor Ir doublets (i,jf =Irl and Ir2). 
Ud^ Um,ab and, Um,dp denote the on-site Coulomb in- 










3 



FIG. 1: Cubic unit cell of the hyperkagome lattice, where red 
sphere stands for iridium atom. The cubic unit cell contains 
12 inequivalent iridium sites. 


teraction within Ir site, the on-site Coulomb interaction 
among three p orbitals within Om oxygen site and inter¬ 
site Coulomb interaction between Ir site and Om oxygen 
site respectively. Based on the strong coupling expan¬ 
sion, the exchange interaction between two Ir J = ^ 
pseudospin doublets is derived from eq. ©• 

By way of Sim,a,aa' and in eq. (|2]), the form of 

the spin interaction between Irl and Ir2 is constrained by 
those additional spatial symmetries applied only for Irl, 

Ir 2 , 01 and 02. A lattice parameter of Na-438i2 dictates 
that these four nearly respect two symmetries; one is a 
bond-centered mirror symmetry (Fig.[2fb)) and the other 
is bond-centered C 2 rotational symmetry (Fig. EJc)). 
The mirror symmetry restricts the Dyzaloshinskii-Moriya 
(DM) vector to be in the mirror plane, while the C 2 ro¬ 
tation constrains the DM vector to be perpendicular to 
the C 2 rotational axis. When the two Ir doublets are 
chosen at Irl and Ir 2 in Fig. [TJ the C 2 rotational axis is 
along ( 0 , 1 , 1 ) and the mirror plane is perpendicular to 
the (0,1, —1), so that the DM interaction vector between 
Irl and Ir 2 is along ( 1 , 0 , 0 ), i.e. D 12 = D 31 = 0 in eq. ([T]). 
Likewise, the symmetric part of the anisotropic exchange 
interaction between Irl and Ir 2 respects G 12 = G31 = 0 
and J 2 = J 3 due to these two approximate symme¬ 
tries. Without any justification, we further assume G23 
in eq. o to be zero. This leads to the following reduced 
spin Hamiltonian for the hyperkagome iridate; 

^(Si,,Y(j + G 0 0 

^ = E M*.- 0 ^ ’(4) 

\ 0 -D J-G \Sj,J 



FIG. 2: (a) A triangle formed by three neighboring iridium 
sites (red spheres; Irl, Ir2 and Ir3 in Fig. [Tj) and their neigh¬ 
boring oxygens (blue spheres; three 01 and one 02). 01 is 
shared by two of three oxygen octahedra whose center ac¬ 
commodate Irl, Ir2, and Ir3 respectively. 02 is shared by 
the three oxygen octahedra. Gubic coordinates (x,y,z coor¬ 
dinates) with which an expression for the d = ^ doublet, 
eq. ©, is defined are depicted by black solid lines with ar¬ 
rows at the iridium sites, (b) Bond-centered mirror operation 
and (c) bond-centered C 2 rotation, under which Irl and Ir2 
are exchanged with each other, while 01 and 02 remain in¬ 
tact. The lattice parameters of the hyperkagome iridate^^ 
suggests that spatial coordinates of 01, 02, Irl, Ir2 and their 
cubic coordinates nearly respect these two symmetries. 


eludes all the effective spin exchange models previously 
derivedThe actual values of J, D and G depend on 
detailed electronic band structure or assumptions. From 
the high-temperature expansion of eq. the Curie- 
Weiss temperature for polycrystalline samples is given 
by 3 J — G. From experimental Curie-Weiss fittings of 
the magnetic susceptibility of powder samplesi^Tii^, we 
focus only on 3 J > G region. 


III. MONTE CARLO SIMULATION 


where ^iij),p) follows the same convention 

as above. This simplified spin model approximately in- 


Monte Carlo (MC) simulation was carried out for a 
12 X spin cluster of the hyperkagome lattice with pe¬ 
riodic boundary condition (L = 6, • • • ,10). We used 
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FIG. 3: Low-temperature magnetic phase diagram in3J — 
G > 0 (T = 0.03). Z 2 phase (red cross points) Ze phases 
(blue double cross points), ferrimagnetic phase (green re¬ 
gion with black filled circle points), Neel phase (grey trian¬ 
gle points), and incommensurate (IC) magnetic phase (white 
region). Three open circles at D = 0 (labeled as (1,2,3)) 
denote high symmetric parameter points; (1) a SU(2) point; 
{J,D,G) = (0,0,—1), (2) AF Heisenberg point; (1,0,0), and 
(3) AF Kitaev point; Four finite-T phase dia¬ 

grams in Fig. |4] (a,b,c,d) are along four arrows with black 
broken lines labelled as (a,b,c,d) in Fig. [3] respectively. 


the Metropolis method combined with an over-relaxation 
(microcanonical) update, where 1 MC step comprises of 
12 X number of single-spin flip trials followed by one or 
two (non-random) sequential applications of microcanon¬ 
ical updato^iiS^i^. To increase an acceptance rate of the 
single-spin flip trial, we impose maximum variation of 
spin to be smaller for lower temperature, \AS\ < T (the 
temperature unit is the square root of -\-D^ The 

microcanonical update comprises of a 7r-rotation of single 
spin around local exchange field created by its four neigh¬ 
boring spins. We apply the 7r-rotation on all the spins in 
sequence. After 30000 ^ 40000 MC steps for the ther¬ 
mal equilibration, physical quantities are measured once 
per every 5 MC steps. Physical quantities are typically 
averaged over 20000 ^ 40000 samples. 

Fig. [3] and Fig. 0] show low-T and flnite-T magnetic 
phase diagrams obtained from the simulation. The low-T 
diagram comprises of Z 2 , Ze, Zg non-coplanar antiferro¬ 
magnetic phases and non-coplanar ferrimagnetic phase, 
all of which are associated with k = 0 ordering, Neel 
phase with k = (0,7r,7r) and incommensurate magnetic 
(IC) phase with k = (a, a, a) with 0 < a < tt. The 
MC simulation does not see any magnetic orderings up 
to T = 0.03 at the two high symmetry parameter points; 
antiferromagnetic (AF) Kitaev point (J = G>0,T = 0) 
and isotropic AF Heisenberg point (J > 0,G = D = 0); 
see appendix A#21i2i 

Z 2 and Ze SB magnetic phases and their stabilities 



FIG. 4: Finite-T phase diagrams along the broken lines in 
Fig. [3l Gritical temperatures (Tc) for Z 2 , Ze SB phase and 
IG phase are all identified with specific heat peak tempera¬ 
tures (black filled circle points). The crossover temperature 
(T*) within the Ze phase is depicted by black upper triangle 
points. Phase boundary between Ze and Z 2 phases (depicted 
by bold purple lines) is of the first order with Zg-like symme¬ 
try (Appendix G). 

can be captured by the Luttinger-Tisza (LT) analy- 
sis j3i^32 rj^Yie analysis begins with the Fourier series of 
the quadratic spin Hamiltonian; 

H = (5) 

with jf, m sublattice index {j^m = 1, • • • ,12) and (3 
spin index. Spin ordering in a magnetic ground state is 
specified by a lowest eigenmode of the 36 by 36 Hermitian 
matrix H{k)^ provided that the eigenmode is real-valued 
and satisfies a fixed norm condition. The fixed norm 
condition requires that the norm over spin index is same 
for different sublattices. A magnetic structure of the Z 2 
phase is given by the lowest eigenmode of H{k) at /c = 0, 
in which spin moment is ordered transverse to the G 2 
rotational axis at respective site. The Z 2 magnetic phase 
is essentially same as the so-called ‘canted windmil phase’ 
discussed in a previous worb^. 

The Ze SB magnetic phase is characterized by the dou¬ 
bly degenerate eigenmodes at k = 0 (|0i) and 102 )), 
which form the 2D irreducible representation oi H(k = 
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FIG. 5: Distribution of (mi,m 2 ) measured at different tem¬ 
perature and at (J, D, G) = (0.50, —0.22, —0.84) with L = 9. 
(a) T > n {T = 0.8 (pink), 0.64 (red), 0.52 (blue), 0.36 
(grey) from center to outside while Tc = 0.691). (b) T < 

(T = 0.24 (red), 0.04 (black) from center to outside). Six 
dotted lines with the Ze symmetry in (b) specify the easy- 
axis directions within the 2D plane. We choose (J, D, G) = 
(0.50,-0.22,-0.84), because this point is proximate to the 
possible candidate parameter point for Na-438^. 


0). Any linear combination of these two cannot satisfy 
the fixed norm condition; moreover, they break the norm 
condition in the U(l) symmetric way. Meanwhile the 
condition on averaged spin moment can be relaxed at fi¬ 
nite temperature due to the thermal fluctuation. In fact, 
the doubly degenerate lowest eigenmodes of H{k) ap¬ 
pear at /c = 0 in the Ze phase region, and spins become 
condensed into a plane subtended by these two below a 
critical temperature Tc (determined by the specific heat 
peak). Fig. [5] shows a distribution of a projection of the 
12 spins onto the 2D plane, i.e. (mi,m 2 ) defined by 

- y (6) 

j,Oi 

where j is the site index and a is the spin index. For 
T > Tc, the projection is accumulated at origin ((^) = 0 
with ^ = mi + im 2 ). For T < T^ the amplitude devel¬ 
ops continuously ((|^|) 7 ^ 0). The simulation on a finite- 
size system suggests another characteristic temperature 
T^{< Tc), above which the distribution has the U(l) sym¬ 
metry {{0) indefinite with T = |T|e^^) but below which 
it acquires an additional Zq structure ((cos 6 ^) 7 ^ 0 with 
T = |T|e^^). 


IV. CLASSICAL ORDER BY DISORDER 
MECHANISM 


The Ze anisotropy at the lower temperature (T < T*) 
can be understood from the classical order by disorder 
mechanism. To see this, let us expand H{k) near k = 0 
with respect to small k and derive an effective Hamilto¬ 
nian in the basis of |0i) and | 02 ) using the k p perturba¬ 
tion. On the second order in k, the effective Hamiltonian 


in the 2D plane subtended by |0i) and \(j) 2 ) is ; 

TJ-2X2/1 \ ^ I t>2x2/i \ I ton{k) ■ tno{k) 

W = eo + {k) + 2^ ----- 


eo is the lowest eigenenergy of H{k = 0) to which |0i) 
and \(j) 2 ) belong to. is a higher eigenenergy of H{k = 
0). is a 2 by 2 block of AH{k) = H{k) - 

H{k = 0) in the basis of \(pi) and 102 ); = 

{(l)i\AH{k)\(l)j). ton{k) is a 2 by V block of AH{k) 
which connects \(j)i) and | 02 ) with the V-fold degenerate 
higher energy eigenstates of H{k = 0) belonging to 

(m = l,--- ,Ar); [ton(fe)]i™ = (0i|Aif(fe)|0);))). 

A form of the 2 by 2 effective Hamiltonian thus ob¬ 
tained is constrained by the point group symmetry. To 
see this, note first that, in the 2D irreducible representa¬ 
tion of II{k = 0 ), ^ 2 ®! •C' 2@55 C' 2@2 •C' 2@4 Uud C 2@3 •C^ 2@6 
are the identity operation; choosing any direction within 
the 2D plane will not break any of these point group 
symmetries. These three require that the Hamiltonian 
is quadratic in small k. Namely, these symmetry oper¬ 
ations change the sign of k, C 2 @i • C 2@5 {kx,ky,kz) = 

{kx, ky, kz), C^ 2@2 ' C^2@4 {kx, ky, kz) — ( kx, ky, kz) 

and C 2@3 • C2m {kx, ky, kz) = {-kx, ky, -kz), while they 
are identity operation in the 2 D plane. Besides, iT*(/c) = 
H{-k) SO that the effective Hamiltonian is real-valued 
within the second order in k. It takes the form of 

H^g'^{k) = d3{k)crs + di{k)cri + do{k)ao + 0{k^) ( 7 ) 

in the basis of \(j)i) and | 02 ); ^ 3 | 0 i/ 2 ) = i| 0 i/ 2 )- Here 
real-valued dj{k) is quadratic in k. 

On the other hand, C 2 @i, € 2^2 and C 2 m are mirror 
operations in the 2D plane, where respective mirror axes 
form 60 degree angle against one another. Let us call a 
mirror axis of C 2 @i to be along 102 ) and call its orthog¬ 
onal vector as |0i). In this basis, three operations are 
represented by; 



These mirror symmetries constrain the form of dj{k) in 
eq. ([7]); 


C 2 @i ■ Hlg‘^{kx,ky,kz) ■ C 2&1 = H^g‘^{-k:c,kz,ky), 

Cy2@2 * -^eff ky^ kz) * C2@2 — -^eff ^ 3:5 kz)-, 

C2@3 ■ Hl^\k,,ky, kz) ■ C2@3 = H^s\kz, -ky, k,). 

These three conditions give out 

Hle‘^{k) =k‘^\d{-2kl + kl + kl)(T3 

- \/^d{ky-k‘l)ai ++ 0{k^), (8) 
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FIG. 6: The lowest magnetic energy of finite-momentum spin 
configuration around \0) = cosO |0i) + sin^ | 02 ) is shown as 
a function of ^ for d > 0. The vertical axis denotes Ee^min = 
minf^Ee^k with Ee^k = {k)\0). 


is an antiferromagnetic (AF) phase without any finite 
off-diagonal elements of magnetic susceptibility tensor; 
Xxy = Xyz = = 0. The last one in the set, C 2 @i • T, 

connects y and z while it connects x with neither y nor 
2 ;; the phase shows spin anisotropy in diagonal suscepti¬ 
bility, 

Xxx 7^ Xyy — Xzz- (9) 

When 1 ^ 2 ) (or its 5 counterparts; 0 = is 

chosen, the configuration breaks C 2 m and C 2@3 , but is 
invariant under C 2 @i; {C 2 m * <^ 205,^^202 * C 2 m,C 2 m * 
^^ 206 : ^^ 2 ®!}- Accordingly, the phase is an AF phase 
without off-diagonal susceptibility tensor and with Xxx 7 ^ 
Xyy = the presence of large lattice-spin coupling, 

both of these two Zq SB phases are accompanied with 
an uniaxial lattice distortion differentiating x coordinate 
from y and z coordinates. 


with k = \k\ and k = kjk. Here d and A{> 2\d\) depend 
on a microscopic parameter. They satisfy A — 2|d| > 0 
in the Zq phase region, where the lowest eigenmodes of 
H{k) are at /c = 0. 

A magnetic energy for a linear combination of |0i) 
and 1 ^ 2 ) with finite (but small) momenta can be cal- 
culated from as Ee^k = {e\H^^^{k)\e) with 

|6*) = cos0|</)i)+sin6*|^2)- The energy for such a spin con- 
figuration can be further minimized with respect to k . 
The minimized energy has 6 minima at^ = 0,±^,±^,7r 
for d > 0, while at 0 = ±-|, for d < 0 (Fig. [6]). 


The optimal directions of k for the former 6 minima are 
along kz^ kx{ky)^ ky{kx)^ kz axes respectively. Due to the 
entropy effect, either one of these two types of six minima 
is selected as the ordering direction at finite temperature. 

To see this entropy effect, one can also calculate the 
free energy which includes spatial fluctuation around \0) 
within the second order in small k: 


m = 


7^4 


E' 


,—l3E0^k 


= — In 
2/3 




with Bx = A — 2dcos20, By = A — 2dcos{20 — and 
Bz = A — 2dcos{20 + ^). At T 7 ^ 0, F{0) has indeed 
six minima at ^ = 0 , ±^, 7 r for d > 0 , while at ^ = 

ford< 0 . Consistently, the MC simulation 
found that the system at the lower temperature prefers 
either one of these two types of six minima (Fig. Efb)). 

When |(/)i) (or its 5 counterparts; 0 = 0,±^,±^) is 
selected by the classical order by disorder, the spin con¬ 
figuration breaks all C 2 @j {j = 1,2,3) but is invariant un¬ 
der C 2 @i -T; the generator of the reduced magnetic point 
group is {C 2 @l • C 2 @ 5 , C 2@2 • C2@4, C 2@3 ' ^ 206 , C 2 @l 'T}. 
The first three in the set change the sign of spin mo¬ 
ments, i.e. C 2 @i ' C 2@5 {Sx-)Sy^Sz^ — {Sx-) Sy^ *^ 2 :), 
C2@2 ' C2@4 {Sx^Sy^Sz^ — ( Sx^ Syi^z) S'Hb C2@3 
C 2 m {Sx,Sy,Sz) = {-Sx,Sy,-Sz). Thus, the phase 


V. EMERGENT U{1) SYMMETRY AROUND 
THE CRITICAL POINT 

The U(l) to Zq crossover behavior around T* in the Zq 
SB magnetic phase can be crudely captured by the six 
states ferromagnetic Potts model in the 3D lattice. Low- 
T physics of the Potts model are controlled by 3D XY 
fixed point, T = 0 Nambu-Goldstone (NG) fixed point 
and T = 0 fixed point with large Zq term j^^i^^ The criti¬ 
cal point (T = Tc) belongs to the XY fixed point, below 
which all the renormalization group (RG) flow goes to the 
large Zq fixed pointThe U(l) to Zq crossover below Tc 
stems from the dangerously irrelevant behavior of the Zq 
anisotropy term j^^i^^ The anisotropy term is renormal¬ 
ized to smaller value around the XY fixed point, while 
it blows up into larger value around the NG fixed point. 
Due to this behavior, a smaller system near Tc behaves 
as if it has much reduced Zq anisotropy, while a larger 
system far from Tc behaves as a system with enhanced 
anisotropy. 

Fig. [7| shows a set of temperature dependencies of the 
XY order parameter and Zq order parameter; 

2 2,2 

m = nil T m2 

niQ = [(mi + im2)^ + (mi - im2)^]/(2m^), 

calculated for different system sizes (I/ = 5,6,7,8,9,(10)) 
at two different parameter points in the Zq SB 
phase region; = (0.50,-0.22,-0.84) and 

(0.69,-0.53,-0.50). The T-dependences of m and me 
clearly show that, in finite system sizes, the XY order 
parameter and the Zq order parameter start to take fi¬ 
nite values at different temperatures. Note also that the 
former parameter point is proximate to the SU(2) sym¬ 
metry point (J, D, G) = (0, 0,-1); the SU(2) point is ap¬ 
proximately an effective spin model of Na-438, when the 
atomic spin-orbit interaction is larger than the non-cubic 
crystal field splitting energy and when the exchange path 
is mainly mediated by the oxygen ions^. Meanwhile, the 
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0.0 0.1 0.2 0.3 0.4 


FIG. 7: Temperature dependence of XY order parameter 
m and Zq SB order parameter (me) at different parameter 
points; (a) = (0.50,-0.22,-0.84) (b) (J,D,G) = 

(0.69,-0.53,-0.50). 


latter parameter point is far from any high symmetric 
parameter points, playing the role of a good reference 
parameter point (see below). 

The crossover system size and temperature can be eval¬ 
uated from the FNS analysis on the XY order parame¬ 
ter m and Zq order parameter tuq. The scaling argu¬ 
ment^-— suggests that these two follow single-parameter 
scalings; 

m = (10) 

me = L~^ fe{tL^^’'^), ( 11 ) 

with G — and t = |T — Td/Tc. [3 and u denote the 
critical exponents of the 3D XY universality class. Fig. [8] 
and [9] show respective one-parameter scaling fittings at 
the two parameter points. To obtain them, we fixed the 
critical exponents [3 and y to be those for the 3D XY uni¬ 
versality class (z/ = 0.672 and p = 0.348)^, while fine- 
tuned the critical temperature Tc such that all the data 
points for m fall into one scaling function f{x). With 
Tc thus obtained, we further fine-tuned the crossover 
exponent for me, z^e, such that all the data points for 



me*L° 



FIG. 8: Single-parameter scalings for XY order parameter 
m (a) and Zq SB order parameter ttiq (b) at (J, D,G) = 
(0.50,-0.22,-0.84). 


me fall into the one-parameter scaling function f6{x). 
The two fittings demonstrate that all the numerical data 
points for m and me taken from different system size 
(L = 5,6,''' ,10) fall into respective one-parameter scal¬ 
ing functions. The optimal uq is evaluated at two differ¬ 
ent parameter points, = (0.69,-0.53,-0.50), 

(0.50,-0.22,-0.84), as 1.45 ±0.05 and 1.85 ±0.05 re¬ 
spectively. The crossover system size A* and temperature 
T* can be obtained by equating the argument of fei''') 
with unit; = 1 or GA^/^® = 1. With ^ ^ we 

obtain A^ ^ ^^j^d AT^ ^ respectively. 

The former value of the crossover exponent uq (1.45 
±0.05) is consistent with previous estimation in the Potts 
model^i^, while the latter value (1.85 ±0.05) is rela¬ 
tively larger. The discrepancy stems from the presence of 
a high symmetric point at (J, D, G) = (0, 0, —1) near the 
latter parameter point. The symmetric point has a global 
SU(2) symmetry toward which the Zq anisotropy dimin¬ 
ishes. Such a symmetric point plays the role of another 
crossover fixed point for RG flow, changing the crossover 
exponent from that of the simple Zq Potts model. 
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FIG. 9: Single-parameter scalings for the XY order param¬ 
eter m (a) and Zq order parameter me (b) at (J, D,G) = 
(0.69,-0.53, -0.50). 


VI. EFFECTS OF QUANTUM FLUCTUATION 

In the quantum J = 1/2 pseudo-spin case, the tech¬ 
niques used in this paper may not be able to capture 
physics at or near highly frustrated parameter points, 
where all magnetic ordering temperatures reduce to zero 
due to strong spin frustration. Such frustrated points in¬ 
clude the AF Heisenberg point (J, D,G) = (1,0,0) and 
Kitaev point (J, I), G) = hand, a 

comparison between preceding classical calculations^i^ 
and quantum calculations^ii^ on the honeycomb iridate 
spin model suggests that, even for the quantum J = 1/2 
case, physics discussed in this paper may hold true in 
those parameter regions with higher magnetic ordering 
temperatures. Such parameter regions cover most of the 
physically relevant parameter regions (3J — G > 0), in¬ 
cluding one of the ‘candidate’ parameter points, i.e. the 
SU(2) point (J, D,G) = (0,0,—1). In these regions, the 
quantum fluctuation changes the results in a quantitative 
level as discussed below. 

Firstly, the Z 2 -Z 6 phase boundaries in Fig. [3] at T = 0 
are expected to move into the Z 2 phase side due to the 
quantum fluctuation. This is because the flnite-T first or¬ 


der phase boundary between Z 2 and Zq phases is slanted 
in a way that the Zq phases are more stabilized against 
the Z 2 phase by the thermal effect (Fig.[Hb,c)). With the 
Clausius-Clapeyron relation, this indicates that, around 
the phase boundary, the Zq SB states have more low en¬ 
ergy classical spin configurations ‘proximate’ to them¬ 
selves than the competing Z 2 SB state has. Here we 
define a given classical spin configuration “x” to be more 
‘proximate’ to than to “z” when x • y > x • where 
denotes the inner product between two spin configu¬ 
rations with respect to spin index and site index. 

Crudely speaking, the quantum zero-point energy^ 
can be regarded as the second-order energy correction 
to the classical ground state, which comes from the vir¬ 
tual hopping processes between the classical ground state 
|0) and another classical spin states with higher classical 
energies |n); 


^sw 


^ Ko|g»P 

„ En- Eq 




( 12 ) 


Here H' stands for the quantum fluctuation term, com¬ 
prising ‘off-diagonal’ terms such as a|a], a^aj, and a\aj 

with i while at being the Holstein-Primakoff boson 
(creation) operator at the i-th site. The zero-th order 
part Hq comprises ‘diagonal’ terms such as ata^, giv¬ 
ing the classical energies to each classical spin states e.g. 
HQ\n) = En\n). 

The indication of the flnite-T phase boundary between 
Z 2 and Zq phases in combination with eq. uni sug- 
gests that, near the Z 2 -Z 6 phase boundaries at T = 0 
in Fig. Ill^b,c), the quantum zero-point energy of the Zq 
states will be larger than that of the Z 2 state. Namely, 
the Zq states have possibly more low-energy classical spin 
configurations (|n)) connecting with Zq states by the lo¬ 
cal perturbation H' than the Z 2 state does. As a result, 
the Z 2 -Z 6 classical phase boundaries in Fig. |ll^b,c) are 
expected to move into the Z 2 phase side, when the quan¬ 
tum fluctuation is included perturbatively. In fact, one 
can find in a recent literature^ a model calculation on 
a different spin system, whose observation agrees with 
the thoughts given above. One may also notice from 
Fig. [IKa) that the finite-T phase boundary between Z 2 
and Zq phases is constrained in the T = 0 plane. This ex¬ 
ceptional feature is, however, due to the additional sym¬ 
metry at T = 0 (Appendix A). 

Another possible quantitative change could be found 
in the finite-size crossover phenomena in the Zq phase, 
which may be effectively characterized by the (d + 1)- 
dimensional Zq ferromagnetic Potts model. In the effec¬ 
tive model, the quantum effect is taken into account as 
an addition of the imaginary time dimension (+1) to the 
spatial dimension (d = 3). The 4D Zq Potts model ex¬ 
hibits the same kind of finite-size crossover phenomena 
with different crossover exponent^. Moreover, in the 
(3+l)D model, finite temperature leads to a non-trivial 
‘finite-size’ effect along the imaginary time direction, in 
the same way as the finite system size does along the 






9 


spatial direction. Thus, one may even expect that the Ze 
ordered phase accommodates two distinct crossover tem¬ 
peratures below a finite critical ordering temperature Tc; 
one is associated with the spatial fluctuation of spins and 
the other with their temporal fluctuation. 


VII. DISCUSSION 

In this paper, we obtained a comprehensive classical 
magnetic phase diagram for the hyperkagome iridate. We 
clarified the origin of the Ze anisotropy and finite-T or¬ 
dering nature of the Ze phases. Our finite-T classical 
phase diagram suggests that the Ze phases could be fur¬ 
ther stabilized against the competing Z 2 phase by the 
quantum order by disorder. 

Based on the finite-T crossover behavior in the Ze SB 
phase, let us finally introduce a possible phenomenology 
of Na-438 powder samples. Firstly we assume that poly¬ 
crystalline grain size is as small as 1/im ^ 0.1/im and we 
regard that a broad peak in specific heat observed around 
T 35i^ in experiments^^-— corresponds to the onset 
temperature of the Ze SB phase (Tc). The FNS argument 
above claims that, below but near Tc, a system smaller 
than the crossover system size behaves as if 

it has no Zq anisotropy term. Spins in such a small grain 
collectively develop a finite amplitude of <F = be¬ 

low Tc, while the phase of 4> still strongly fluctuates and 
so does that of individual spins. Thus, any local spin mo¬ 
ment seen by probe spins can be averaged to be zero in 
the intermediate temperature regime T* < T < Tc. For 
^ = 3, 5, 7,10 cubic unit cells, the crossover system size 
is evaluated to be A* = 20, 90, 230, 630 cubic unit cells 
for z/6 = 1.85 (0.017, 0.08, 0.20, 0.56 /im with a=8.95 A). 
On further lowering temperature, the correlation length 
^ becomes shorter and so does A*. When A* exceeds 
the grain size on lowering temperature (T < the Zq 
anisotropy becomes prominent and the phase of 4> starts 
to be locked into the six minima. As explained above, 
this locking breaks the point group symmetry, giving rise 
to the uniaxial lattice distortion via spin-lattice coupling. 
Since neighboring grains in polycrystalline sample are ex¬ 
pected to be randomly oriented against one another, the 
locking or its onset around T* with the uniaxial distortion 
will conflict with the grain structure. This may result in 
a ‘configurationally degenerate phase with fluctuating or¬ 
der’ as suggested in the experiments^. 

The author would like to thank Gang Chen, Kenji 
Harada and Xuerong Liu for fruitful discussions. This 
work was financially supported by NBRP of China 
(2015CB921104). 


Appendix A: Finite-T phase diagram at T = 0 

Finite-T phase diagrams at T = 0 are shown in Fig.flQl 
A group of Hamiltonians at T = 0 is symmetric un¬ 
der a ‘Klein’ transformation G^; Hq = G ’ • G 



FIG. 10: Finite-T phase diagrams at T = 0 and 3J — C > 0. 
(a) ^-\-(p<0<7r-\-(p{h)(p<0< with (J, G) = 

(sin cos and ianip = 1. The critical temperatures of the 
Zs AF phase and ferrimagnetic phase are determined from 
specific heat peak (black filled circle points and black dou¬ 
ble crossed points). The crossover temperature within the 
Zs AF phase is determined from temperature dependence of 
sublattice magnetization (black upper triangle points). 


with (J, G) = (sin6>, cos6>). This connects antiferro¬ 
magnetic (AF) and ferromagnetic (F) Heisenberg point 
(J, G) = (±1,0) with two other symmetric points with 
global SU(2) symmetries (T, G) = (0,±1), while leaves 
intact AF Kitaev point (J, G) = F Ki- 

taev point {J^G) = Consistently, the finite-T 

phase diagram at T = 0 is symmetric under this trans¬ 
formation, where the ferrimagnetic phase (f < ^ < f) 
is transformed by G to the Zg antiferromagnetic phase 
(0 < 6> < ^), the ferromagnetic phase (—^ < 6> < 0) to 
the Zg antiferromagnetic phase (f < ^ < x)- "Fhe or¬ 
dering temperatures of respective phases vanish at two 
AF SU(2) points (0 = 0, |^) and two Kitaev points 
{0 = f,^). On the one hand, two F SU(2) points 
{0 = TT,^) are unfrustrated points, where magnetic or¬ 
derings show the maximum transition temperature. No 
0(3) to Zg crossover is observed within 0 < ^ < while 
,for ^ < 0 < 27r, 0(3) to Zg crossover are observed be¬ 
low ordering temperatures for a finite-size system. The 
crossover region becomes widest at the unfrustrated F 
SU(2) point. 


Appendix B: finite size scaling analysis in Zg SB 
phase 

The MC simulation shows that the phase boundary 
between the Z 2 and Zq SB phase is always of the first or¬ 
der (Appendix C), indicating the existence of a bicritical 
point at finite temperature. Below the bicritical point, 
the system exhibits Zg SB magnetic phase with 0(3)- 
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FIG. 11: Temperature dependence of 0(3) order parame¬ 
ter m and Zg SB order parameter mg at different parame¬ 
ter points; (a) {J,D,G) = (0.11,0.0,-0.99) (b) {J,D,G) = 
(0.31,0.0,-0.95). 


type spin fluctuation. This situation can be clearly seen 
at the phase boundary of = 0 by the LT analysis. In 
the phase boundary between Z 2 and Zq phases, H{k) has 
triple degeneracy in its lowest eigenmodes at = 0; two 
are the doubly degenerate lowest eigenmodes from the Zq 
phase region {\(j)i) and | 02 )) and one is the lowest eigen¬ 
mode from the Z 2 phase region (l^s)). These three form a 
0(3) sphere in the 36 dimensional space. On the sphere, 
the fixed norm condition is satisfied only along 8 high- 
symmetry directions; ±(1,1,1), ±(1,1, —1), ±(1, —1,1), 
±(-l, 1,1) ((1,1,1) means |0i) + | 02 ) + l^s)), while oth¬ 
erwise not in general. These directions become strong 
easy-axis directions, giving rise to the cubic anisotropy 
in the 0(3) vector model. 

The ordering nature of the 3D 0{n) vector model with 
the cubic anisotropy is determined by either 3D Heisen¬ 
berg fixed point {n < ric) or cubic fixed point {n > ric). 
Preceding studies including the fifth order 4-e expansion 
conclude that ric is smaller than 3 {ric = 2.89)Being 
consistent with this, the MC simulation on G > 0 and 
D = 0 region does not observe any 0(3) to Zg crossover 
behavior (see Fig. [TOT b)): Below a critical temperature 
(determined by the specific heat peak), the Zg anisotropy 



FIG. 12: Single-parameter scalings for the Zg SB order param¬ 
eter at D = 0 with different parameter points; (a) (J, G) = 
(0.11,-0.99) (b) (J,G) = (0.31,-0.95). 


always becomes prominent in the 0(3) sphere even for 
the smallest simulated system size (L = 6). 

On the one hand, we also observe an unexpected 0(3) 
to Zg crossover region below the bicritical point at G < 0 
and D = 0 (Fig. [TOT a)). Fig [TT] shows a set of tempera¬ 
ture dependences of the 0(3) order parameter m and Zg 
order parameters ruc,^ (/i = 1,2,3) calculated for differ¬ 
ent system sizes and at two different parameter points in 
the Zg SB phase region at G < 0 and D = 0; 


m = y nn\ + 

_ _ {rriy + irUpY + {m^ - irUpY 

Three Zg SB order parameters (/i = 1, 2, 3) show the same 
temperature dependence. Fig. [12] shows one-parameter 
scaling forms for the Zg SB order parameters; ruc^p = 
L~^gc{tL^^^Y' To obtain them, we fixed u and {3 to 
be those of 3D cubic fixed point {u = 0.704 and p = 
0.362)^, while we fine-tuned Tc and Vc to fit data points 
for m and rric^p {jJi = 1,2,3) into g{x) (not shown) and 
gc{x) (Fig. [T2|) respectively. The optimal Vc thus deter¬ 
mined is 1.55±0.05 at the two different parameter points. 
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FIG. 13: Temperature dependencies of Z 2 SB order param¬ 
eter (red cross points), XY order parameter (black filled 
circle points) and specific heat (blue double cross points) 
around the phase boundary between Z 2 SB phase and Ze SB 
phase; (a) [j,D,G) = (0.05,-0.36,-0.93) (b) iJ,D,G) = 
(-0.05,0.25,-0.97) 


The 0(3) to Zs crossover at G < 0 and D = 0 is ap¬ 
parently counterintuitive from the viewpoint of the pre¬ 
vious studies on the 0{n) vector model with the cubic 
anisotropy^^ A simple explanation for this is that the 
cubic anisotropy at the microscopic Hamiltonian level is 
much smaller than Tc and because the largest simulated 
system size {L = 10) is still too small that renormal¬ 
ized (thus enhanced) cubic term is tiny compared to Tc. 
In fact, the cubic term vanishes completely at the un¬ 
frustrated SU(2) point; (J, T,G) = (0,0,—1), while Tc 
takes the largest value at this SU(2) point (see 0 = n 
in Fig. [ToTa)). Besides, a positive scaling dimension of 
the cubic term around the 3D Heisenberg fixed point was 
suggested to be very small.— 


Appendix C: MC evidence for the first order phase 
transition between Z 2 and Ze SB phases 

Fig. [13] shows a set of temperature dependence of the 
Z 2 SB order parameter and that of XY order parameter 


(order parameter for the Zq SB phase) around finite-T 
phase boundaries between Z 2 phase and Zq phase. The 
results show discontinuous changes of these two order 
parameters at the phase boundaries, indicating that the 
transition is of the first order. 

To obtain them, we defined these two order parame¬ 
ters as follows. According to the group theory analysis^, 
H{k = 0) has two ID irreducible representations which 
have the same index table and which correspond to the 
lowest eigenstate of H{k = 0) in the Z 2 SB phase region. 
Let us call their bases as \4>^^) and respectively. 

Note also that H{k = 0) has three 2D irreducible rep¬ 
resentations which have the same index table and which 
correspond to the doubly degenerate lowest eigenstates 
of H{k = 0) in the Zq SB phase region. Call respective 

three sets of doubly degenerate bases as \4>i^) and \ 4 > 2 ^) 
(/ = 1,2,3). 

The order parameter for the Z 2 SB phase is given by 
a linear combination of the following two quantities; 

tm = 


with m = 1,2 and j sublattice index, a spin index. In 
Fig. HSl the temperature dependence of either ti or ^2 
is shown as the Z 2 SB order parameter. The XY order 
parameter for the Zq SB phase is given by a combination 
of the following three quantities; 


with / = 1, 2, 3 and 


with m = 1,2. In Fig. [131 ^Fe temperature dependence 
of one of and is shown as the XY order 

parameter for the Zq SB phase. 
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